eHSP90α in front-line therapy in EGFR exon 19 deletion and 21 Leu858Arg mutations in advanced lung adenocarcinoma

Purpose Extracellular heat shock protein 90 AA1(eHSP90α) is intricately linked to tumor progression and prognosis. This study aimed to investigate the difference in the value of eHSP90α in post-treatment response assessment and prognosis prediction between exon 19 deletion(19DEL) and exon 21 Leu858Arg(L858R) mutation types in lung adenocarcinoma(LUAD). Methods We analyzed the relationship between the expression of eHSP90α and clinicopathological features in 89 patients with L858R mutation and 196 patients with 19DEL mutation in LUAD. The Kaplan-Meier survival curve was used to determine their respective cut-off values and analyze the relationship between eHSP90α expression and the survival time of the two mutation types. The area under the curve (AUC) was used to evaluate the diagnostic performance of biomarkers. Then, the prognostic model was developed using the univariate-Cox multivariate-Cox and LASSO-multivariate logistic methods. Results In LUAD patients, eHSP90α was positively correlated with carcinoembryonic antigen(CEA), carbohydrate antigen 125(CA125), and carbohydrate antigen 153(CA153). The truncated values of eHSP90α in L858R and 19DEL patients were 44.5 ng/mL and 40.8 ng/mL, respectively. Among L858R patients, eHSP90α had the best diagnostic performance (AUC = 0.765), and higher eHSP90α and T helper cells(Th cells) expression were significantly related to shorter overall survival(OS) and worse treatment response. Also, high eHSP90a expression and short progression-free survival(PFS) were significantly correlated. Among 19DEL patients, CEA had the best diagnostic efficacy (AUC = 0.734), and CEA and Th cells were independent prognostic factors that predicted shorter OS. Furthermore, high CA125 was significantly associated with short PFS and poor curative effect. Conclusions eHSP90α has a better prognostic value in LUAD L858R patients than 19DEL, which provides a new idea for clinical diagnosis and treatment. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-024-12573-3.


Introduction
Lung cancer has been shown to have the highest cancerrelated mortality rate of all cancers [1].Furthermore, the proportion of lung adenocarcinoma (LUAD) is the largest and still increasing [2].Lung cancer treatment mainly includes surgical, chemotherapy, radiotherapy, and immunotherapy [3].With the significant increase in lung cancer pathogenesis research and precision therapy, molecular targeted therapy based on lung cancer driver mutations has dramatically enhanced the overall survival(OS) rate of patients [4].Epidermal growth factor receptor (EGFR) mutations occur mainly in LUAD.EGFR tyrosine kinase inhibitor (EGFR-TKI) sensitizing mutations are most commonly found in exon 19 deletion (19DEL) or exon 21 Leu858Arg (L858R) [5].Moreover, some studies have found differences in the efficacy of targeted therapy for different EGFR mutation types, with some mutations being more responsive to treatment than others [6,7].This difference is crucial for LUAD patients' precise treatment and prognosis.
The heat shock protein 90 (HSP90) family comprises a highly conserved chaperone protein with a molecular weight of approximately 90 kDa, widely distributed in mammals.HSP90 has different isoforms in various intracellular spaces, with heat shock protein 90 AA1(HSP90α) and heat shock protein 90 AB1(HSP90β) primarily distributed in the cytoplasm.HSP90α is the stress-inducible type, and HSP90β is expressed constitutively [8][9][10].Besides its intracellular localization, HSP90α can also be excreted into the extracellular environment, referred to as extracellular HSP90α (eHSP90α).As a chaperone protein, HSP90α has numerous client proteins and is therefore involved in many critical activities, such as inflammation regulation, apoptosis, and immunity [11].Previous research found that the expression of eHSP90α is higher in cancer patients than in healthy populations, making it a reliable prognostic marker for several cancers [12].For instance, Han et al. demonstrated that eHSP90α could be an effective diagnostic biomarker for liver cancer and predict patient response to surgery [13].Similarly, in lung cancer(LC), Shi et al. found it can be a valid diagnostic biomarker and can expect the response to chemotherapy [14].Besides, Huang et al. found that high expression of eHSP90α is associated with the poor efficacy of chemotherapy and prognosis in Small Cell Lung Cancer(SCLC), and the area under the diagnosis curve of eHSP90α against SCLC appears to be 0.791, with outstanding sensitivity and specificity [15].
However, we have yet to find any studies exploring whether eHSP90 is different under different EGFR mutations in LUAD.This study aims to use eHSP90α to predict the prognosis of patients with 19DEL and L858R mutations in LUAD patients and explore its different prognostic values.In addition, we further combined eHSP90α with other serum markers to develop a predictive nomogram model to assess prognosis accurately.

Patient inclusion and exclusion criteria
This study included 93 patients with L858R in LUAD and 220 patients with 19DEL mutation admitted to the Department of Respiratory Oncology, Guangxi Medical University Cancer Hospital from July 2008 to July 2021.Four L858R patients and 16 19DEL patients were excluded due to incomplete clinical information.Eight 19DEL patients with exon 20 mutations were also excluded (Fig. 1A).
The inclusion criteria for this study were as follows: (1) Pathological diagnosis of LUAD; (2) Obtained blood samples; (3) Confirmed EGFR mutations of 19DEL or L858R only; (4) Availability of tumor markers and immune cells after admission.(5) Response assessments were conducted using chest and abdominal computed tomography scans (CT scans) every two treatment cycles, following the Response Evaluation Criteria in Solid Tumors (RECIST version 1.1).Each patient's response was categorized into one of the following groups: responders, which encompassed cases of complete response (CR), partial response (PR), and stable disease (SD), and non-responders, which included instances of disease progression (PD).Exclusion criteria for this study were: (1) Incomplete clinical information such as tumornode-metastasis (TNM) stage; (2) They also carried other EGFR mutations or have not been tested for EGFR gene mutations; (3) Patients who also had other oncological diseases.;(4) Incomplete follow-up information (Table 1).
Our outcome overall survival(OS) was based on all causes of death, while progression-free survival(PFS) was based on time for people who relapsed.Follow-up for all participants ended in July 2021.Among them, the median follow-up time of L858R mutation patients was 15 months, and for 19DEL mutation patients, it was 29 months.

Information collection
This prospective study collected peripheral blood data of patients diagnosed with LUAD who were detected as 19DEL or L858R mutation, including age, gender, TNM staging, carbohydrate antigen 125 (CA125), carbohydrate antigen 153 (CA153), carbohydrate antigen 199 (CA199), carcinoembryonic antigen (CEA).In addition, immune cells, T lymphocytes (T cell), T helper cells (Th), suppressor T cells (Ts), T helper cells/ suppressor T cells (Th/Ts), natural killer (NK) cells and B lymphocytes (B cell) were included.All blood samples were collected from patients for the first time after admission.

Detection of serum HSP90α
HSP90α was measured using a Yantai Protgen Biotechnology Development Co., Ltd., Shandong, China assay kit.Patient blood samples, collected in EDTA anticoagulant tubes in the early morning, were centrifuged at 3000 rpm for 15 min at 4 °C to remove particles.The assay kit was equilibrated at 37 °C for 30 min before use.Plasma samples were diluted 20 times with the diluent.Each assay plate well received 50µL of standard, plasma sample, and HSP90α marker solution.Incubating the plate at 37 °C for 60 min.After incubation, the wells were washed multiple times with the washing solution.Color developers A and B were added sequentially, and the plate was incubated at 37 °C in the dark for 20 min.The color development was then stopped by adding a stop solution.The optical density (OD) values were measured at 450 nm/620 nm within 10 min of color development termination.The HSP90α content in plasma samples was calculated based on these OD values.

Nomogram, lasso regression, logistic regression, and cox regression
The nomogram consists of a set of parallel lines, each representing a variable and a scale along each line.Lasso Regression is a type of linear regression that uses regularization methods to prevent overfitting.In traditional linear regression, the model tries to fit the data as closely as possible, which can lead to overfitting and poor performance on new data.Logistic regression is defined as: p = \frac {1} {1 + e^{-z}}.where 'p' is the predicted probability, 'z' is the weighted sum of the input features, and 'e' is the mathematical constant of approximately 2.71828.The Cox regression model estimates the hazard ratio, which is the ratio of the hazard rates between two groups while adjusting for other covariates that may affect the outcome.

Statistical analysis
Data analysis was performed using SPSS25.0software.Clinical baseline data were presented as medians and interquartile ranges.Spearman rank correlation analysis was used to analyze the correlation between eHSP90α and other clinical biomarkers.The ROC curve was used to evaluate the diagnostic efficacy of all included clinical indicators for the prognosis of the disease.LASSO and multivariate logistic regression analyses were used to analyze the relationship between response to treatment and eHSP90α.Univariate and multivariate Cox regression analyses analyzed the relationship between eHSP90α and OS and PFS.Peripheral blood indicators less than 0.05 were included in the multivariate analysis.The best cut-off value of each index was calculated using R4.03 software, and the Kaplan-Meier(K-M) survival curve was drawn.Nomogram software was used for validation.P values < 0.05 were considered statistically significant.
Furthermore, we analyzed the K-M survival curves of OS for two mutations and tested them with log rank.The K-M survival curve shows that patients with 19DEL mutation predicted a worse prognosis (Fig. 1D), and the median survival time of L858R mutation patients was 43 months; the median survival time of patients with 19DEL mutation was 63 months.Also, statistics show that the OS of L858R patients was lower than that of 19DEL patients (Fig. 1C).

Correlations of eHSP90α and clinical index
CEA is the most common tumor marker for LUAD, and we analyzed the correlation between eHSP90α and CEA in LUAD and with other clinical biomarkers (Table 3).Through analysis, we found that eHSP90α was positively correlated with CEA, CA125, and CA153 among all LUAD patients we recruited.When differentiating mutations, eHSP90α was only positively associated with CEA in patients with L858R mutation.In patients with 19DEL mutation, eHSP90α was positively correlated with CEA, CA125, and CA153.eHSP90α may have good diagnostic and prognostic value, but there may also be some differences between the two mutations.

Inclusion criteria
Exclusion criteria They have pathologically diagnosed LUAD.
Clinical information such as TNM stage was incomplete.Blood samples could be obtained.
They also carried other EGFR mutations or has not been tested for EGFR gene mutations.19DEL or L858R was the only EGFR mutation they carried.
Patients who also had other oncological diseases.Tumor markers and immune cells were availabile.
Follow-up information was incomplete.Response assessments were conducted using chest and abdominal computed tomography scans (CT scans) every two treatment cycles, following the Response Evaluation Criteria in Solid Tumors (RECIST version 1.1).

Diagnostic value of biomarkers in LUAD
To better study the influence of clinical and biological indicators on the prognosis of LUAD patients, we made Receiver operating characteristic(ROC) curves of each physical index with OS as the outcome.We calculated the area under the curve(AUC) to judge its diagnostic efficiency (Table 4).Through analysis, we found that the top three AUC values in all our included patient data were CEA, CA125, and eHSP90α, with 0.718, 0.668, and 0.644, respectively.In L858R mutant LUAD patients, the top three places calculated by the area under the curve were eHSP90α, CEA, and CA125, and their AUC values were 0.765, 0.682, and 0.63, respectively.In patients with 19DEL mutation, the top 3 areas under the curve were CEA, CA125, and CA153, and their AUC values were 0.734, 0.684, and 0.673, respectively, while eHSP90α had an AUC of 0.591 and a Youden index of 0.209.Suppose a single biomarker is

Value of eHSP90α in clinical prognostic assessment in patients with L858R mutation and 19DEL mutation
Then, we explored the different expressions of eHSP90α in Alive and Dead groups in OS in patients with L858R and 19DEL mutations, respectively.Statistics show that in L858 patients, those who eventually died had higher eHSP90α levels than the non-dead group, whereas in 19del, there was no difference in eHSP90α levels between the two groups (Fig. 1B).
The K-M survival analysis was performed to determine the cutoff points of the index.The cut-off value of eHSP90α for L858R mutation with OS was 44.5ng/mL, and for 19DEL mutation was 40.8ng/mL.Kaplan-Meier survival analysis showed that L858R mutation patients with high eHSP90α levels had poor OS.However, this trend could not be detected in19DEL mutation patients (Fig. 1E).Other clinical biomarkers were also subjected to K-M survival analysis, and the indicators that have predictive value on prognosis are shown in Figs.1E and 2.
We performed COX univariate and multivariate analyses for OS outcomes to determine which factors are independent in the prognosis of patients with L858R and 19DEL mutation (Tables 5 and 6).In multivariate analysis, we found that in patients with L858R mutation, high eHSP90α expression and high Th levels were short OS significantly associated (eHSP90a, P = 0.0305 and Th, P = 0.0282); In patients with 19DEL mutations, high CEA and Th levels had independent indicators of the prognosis of their short OS (CEA, P = 0.0150, and Th, P = 0.0392).Therefore, we included eHSP90α and Th, CEA, and Th as reference factors for two different mutations, constructing their nomograms to predict the OS of two other mutations (Fig. 3).

Indicator value of eHSP90α in PFS and response evaluation in patients with L858R mutation and 19DEL mutation
Compared with OS, PFS and treatment response (RECIST version 1.1) reflect the short-term survival benefits of patients more.Therefore, PFS was first used as the outcome for univariate and multivariate analysis of Cox (Tables 5 and 6).In multivariate analysis, we found a significant relationship between high eHSP90α expression and short PFS in patients with L858R mutation (eHSP90α, P = 0.0173).In patients with 19DEL mutations, high CA125 is an independent indicator of the prognosis of PFS (CA125, P = 0.000575).

Discussion
EGFR is a transmembrane glycoprotein distributed on the surface of epithelial cells and consists of extracellular ligand-binding regions, transmembrane regions, and intracellular tyrosine kinase binding regions.Studies have found that patients with EGFR gene mutations are relatively more likely to have metastasis.The most common is bone metastasis and brain metastasis [16].In recent years, targeted therapy for patients with EGFR mutations has achieved significant clinical efficacy, significantly prolonging survival compared with traditional chemotherapy regimens [17].19DEL and L858R are the two most common sensitive mutations in EGFR-TKI.Data from multiple trials suggest that the efficiency of different treatment strategies may vary depending on various EGFR mutation states, particularly between 19DEL and L858R mutations.For example, studies have found that patients with 19DEL have a higher PFS than patients treated with TKIs than patients with L858R [18].It has also been found that OS in patients with 19DEL is higher than in patients with L858R in our study [7].HSP90α is the only isoform detected in plasma.In addition to participating in various cell life activities such as apoptosis, it has also been used as a clinical biomarker in multiple cancers.Several studies have shown that it has good diagnostic performance and prognostic value as a marker for cancers [19][20][21].Studies have long found that eHSP90a can assist in  diagnosing LC and is associated with its progression [14,22].Still, no one has explored the expression and diagnostic efficacy of eHSP90a in different subtypes of EGFR mutations in LUAD.Our previous study found that in SCLC, the diagnostic efficacy of eHSP90 as a prognostic evaluation and diagnostic marker is better than NSE, a biomarker of classical SCLC [15]; In this study, we found a positive correlation between eHSP90α and CEA.In the evaluation of diagnostic efficacy, we found that among the patients with L858R mutation, eHSP90α had the best diagnostic efficacy (AUC = 0.765) among all the indicators we included in the evaluation and was only slightly inferior to CEA in a single specific index; In the data of patients with 19DEL mutation, we found that the performance of eHSP90α is somewhat insufficient, and the biological indicator with the highest diagnostic efficacy is CEA (AUC = 0.734), despite this, CEA is still not as sensitive as eHSP90α.Another study also found that eHSP90α is a valuable predictor of response to early chemotherapy and positively correlated with tumor remission after chemotherapy in NSCLC [23].This study found that high eHSP90α was closely related to short OS, PFS, and progression after treatment in L858R mutation patients.Among 19DEL mutation patients, high eHSP90α can only predicted poor treatment response.eHSP90α is a highly sensitive pan-cancer marker; its specificity may not be enough, so this paper for two different mutation types, using eHSP90α as the basis, combined with other markers to construct their corresponding prognostic models.This study also found that high Th cells were closely related to poor OS in 19DEL and L858R mutation patients.Previous works of literature found that NSCLC patients in the advanced group were evidently lower in the expression of CD4 + but markedly higher in the expression of CD8 + in peripheral blood than in the early group [24].Another study found that LUAD has higher CD4 + expression than other subtypes in NSCLC [24,25].Moreover, high expression of CD4 + may indicate a better immune response and should indicate better patient survival.Thus, the mechanism still needs to be explored.
In summary, this study found that patients with L858R mutation had worse OS than those with 19DEL mutation.In patients with L858R mutation, eHSP90α is closely related to short OS, PFS, and progression after treatment.As a new biomarker, eHSP90α is highly valued in evaluating its prognosis.However, in patients with 19DEL mutation, the indicative value of eHSP90α is relatively limited, and it only has a specific, meaningful value for the response to treatment.However, this study also had some limitations.First, this is a single-center study; the data may had some selection bias.Second, the accuracy of our nomogram should be evaluated through external validation, which will help to assess whether

Fig. 1
Fig. 1 Workflow diagram and survival curves of eHSP90a expression differences, survival in lung adenocarcinoma (LUAD) EGFR exon 19 Deletion(19DEL) and 21 Leu858Arg(L858R) mutations patients, and Kaplan-Meier Survival curve of eHSP90a and overall survival(OS) in patients with LUAD L858R mutation.(A) Workflow.(B) The scatterplot shows the difference in expression of the continuous variable eHSP90a in the OS outcome state of the 19del and L858R mutations.(C) The scatterplot shows the OS of 19DEL and L858R mutation.(D) Survival of 19DEL and L858R mutation patients.(E) Kaplan-Meier Survival curve of eHSP90a and OS in patients with LUAD L858R mutation

Fig. 3 Fig. 4
Fig. 3 Nomograms based on OS in patients with LUAD 19DEL and L858R mutations.The prognostic nomogram for OS is based on the prognostic scores of eHSP90α and other factors in (A) L858R and (C)19DEL patients.The prognostic nomogram for PFS is based on the prognostic scores of eHSP90a and other factors in (B) L858R and (D) 19DEL patients

Table 2
Baseline on L858R and 19DEL categorical variables in our study used to diagnose LUAD among the L858R mutations.In that case, the highest specificity is CEA, the highest sensitivity is Ts, and the highest Youden index is eHSP90α, with a Youden index of 0.505.In contrast, among the 19DEL mutations, the highest specificity is T cell, the highest sensitivity is CA.153, and the highest Youden index is CEA, with a Youden index of 0.397.So, the CEA index has the best diagnostic value in patients with 19DEL mutation, while eHSP90α has a better diagnostic value in L858 mutation patients in LUAD.

Table 5
Cox regression analysis on the L858R OS and PFS of LUAD

Table 6
Cox regression analysis on the 19DEL OS and PFS of LUAD Th: T helper cells; Ts: suppressor T cells; B cell: B lymphocytes; Th/Ts: T helper cells/ suppressor T cells; NK: natural killer; TK1: thymidine kinase 1; OS: overall survival; PFS: progression-free survival; CI: confidence interval; HR: hazard ratio; OS: overall survival; PFS: progression-free survival